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Abstract 

If the soft Supersymmetry (SUSY) breaking masses and couplings are complex, 
then the associated CP-violating phases can in principle modify the known phe- 
nomenological pattern of the Minimal Supersymmetric Standard Model (MSSM). 
We investigate here their effects on Higgs boson production in the gluon-gluon fusion 
mode at the Tevatron and the Large Hadron Collider (LHC), by taking into account 
all experimental bounds available at present. The by far most stringent ones are 
those derived from the measurements of the Electric Dipole Moments (EDMs) of 
fermions. However, it has recently been suggested that, over a sizable portion of 
the MSSM parameter space, cancellations among the SUSY contributions to the 
EDMs can take place, so that the CP-violating phases can evade those limits. We 
find a strong dependence of the production rates of any neutral Higgs state upon 
the complex masses and couplings over such parts of the MSSM parameter space. 
We show these effects relatively to the ordinary MSSM rates as well as illustrate 
them at absolute cross section level at both colliders. 



1 Introduction and plan 



The soft SUSY breaking parameters of the MSSM can well be complex. Even in the 
absence of flavour non-conservation in the sfermion sector, the Higgsino mass term, the 
gaugino masses, the trilinear couplings and the Higgs soft bilinear mass need not be real. 
Assuming universality of the soft gaugino masses at the Grand Unification (GUT) (or 
Planck) scale, the effects of complex soft masses and couplings in the MSSM Lagrangian 
(see Appendix A) can be parametrised in terms of only two independent phases [|1|, 0^ 
and (pAi associated to the (complex) Higgsino mass term, /x, and to the trilinear scalar 
coupling y4[|, respectively. In other terms. 

Their presence is a potentially dangerous new source of violation of the CP-symmetry in 
the MSSM. But their size can in principle strongly be constrained by the measurements 
of the fermionic EDMs (mainly, of electron and neutron) and several analyses [Q] have 
indicated that 0^ and (pA must be small in general. However, recent investigations p, ^ 
have shown that, in a restricted but still sizable part of the parameter space of the 
MSSM, the bounds drawn from the EDM measurements are rather weak, so that such 
phases can even be close to 7r/2. This is a consequence of cancellations taking place 
among the SUSY loop contributions to the EDMs. Although, in order to be effective, 
these require a certain amount of 'fine-tuning' among the soft masses and couplings 0, 
it has recently been suggested that such cancellations occur naturally in the context of 
Superstring models 0. If the SUSY loop contributions to the EDMs do vanish, then, 
as emphasised by the authors of Ref. SUSY parameters with large imaginary parts 
may have a non-negligible impact on the confrontation of the MSSM with experiments. 
In particular, many of the SUSY (s)particle production and decay processes develop a 
dependence on 0^ and 0^, so that, in view of the importance of searches for New Physics 
at present and future accelerators, their phenomenology needs a thorough re-investigation. 

Various sparticle processes including the effect of such phases have recently been con- 
sidered. For example, neutralino and chargino ||^ production at LEP, at the CERN 
LHC as well as at future electron-positron linear colliders (LCs) [TU, rT|. Direct 



CP-asymmetries in decays of heavy hadrons, such as 5 — * + 7, -B ^ + 7 and 
B — > Xsl^l~ , have been investigated in the context of the Supergravity inspired MSSM 
(M-SUGRA): in Refs. [|12[, 1T|] and |T|], respectively. 



In this paper, we are concerned with the Higgs sector of the MSSM. Here, although 
the tree-level Higgs potential is not affected by the CP- violating phases, since Mfj^ , Mfj^ 
and tan/? (the mass parameters of the two Higgs fields and the ratio of their vacuum 
expectation values (VEVs), respectively) are real and /i enters only through it should 
be noticed that this is no longer the case if one includes radiative corrections. In their 

-'^For simplicity, hereafter (except in the Appendices), we assume A = A„ = at the electroweak 
(EW) scale, i.e., 0{Mz), where u and d refer to all flavours of up- and down-type (s)quarks. 
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presence, one finds |]l5i |T6[| that the three neutral Higgs bosons can mix and that their 
effective couplings to fermions can be rather different at one loop. However, for the MSSM 
parameter space that we will consider here, such corrections turn out to be negligible, of 
the order of just a few percent, as compared to those induced at the lowest order by the 
CP-violating phases in the squark sector. There are some reasons for this. First of all, 
the induced radiative corrections to the Higgs-quark-quark vertices can be parametrised 
in terms of the mass of the charged Higgs boson, Mh± (~ Mjfi). Then, one can verify 
that they essentially depend only upon the input values given to \A\ and MgusY (the 
typical mass scale of the SUSY partners of ordinary matter). Here, we will mainly be 
concerned with trilinear couplings in the range \A\'^ 700 GeV, Higgsino masses of 
the order of 600 GeV or so, and Msusy — 300 GeV. According to the analytic formulae 
of Ref. [0, in the above MSSM regime, one finds negligible corrections to the tree-level 
h^ti and h^hb couplings. (Similarly, for the case of the lightest Higgs boson mass.) In 
contrast, the strength of the H^ti vertex can significantly be modified for not too heavy 
masses of the charged Higgs boson (say, Mh± ~ M^o < 200 GeV) and rather large values 
of |A| and \ix\ (typically, |y4| ~ \^\ ~ 2 TeV), a region of parameter space that we will 
avoid, whereas that of the H%h one is generally small because we shall limit ourselves to 
the interval 2 ^ tan (3 ^ 10. As for papers allowing instead for the presence of non-zero 



values of 0^ and/or (pA and studying the Higgs sector, one can list Refs. |]T3 and fT^. The 
first publications deals with decay rates whereas the second one with Higgs production 
channels probable at a future LC. In such papers though, no systematic treatment of 
the limits imposed by the EDM measurements was addressed. Such effects ought to be 
incorporated in realistic analyses of the MSSM Higgs dynamics. 

Here, we have necessarily done so, since it is our purpose to study the effect of fi- 
nite values of 0^ and/or 0^ on Higgs production via the gg —>■ $° channel |19], where 



$0 = if 0^ j\o represents any of the three neutral Higgs bosons of the MSSM. (A 

preliminary account in this respect was already given by the authors in Ref. [0.) These 
processes proceed through quark (mainly top and bottom: i.e., t and b, see Fig. 1) and 
squark (mainly stop and sbottom: i.e., ii, t2 and 6i, 62, see Fig. 2, each in increasing order 
of mass) loops, in which the (s)fermions couple directly to . Needless to say, as in the 
MSSM the lightest of the Higgs particles is bound to have a mass not much larger than 
that of the Z boson, M^, much of the experimental effort at both the Tevatron (Run 2) 
and the LHC will be focused on finding this Higgs state, h^. In this respect, we remind 
the reader that direct Higgs production via gluon-gluon fusion is the dominant mechanism 
over a large portion of the MSSM parameter space at the LHC and a sizable one at the 
Tevatron 

The plan of the paper is as follows. The next Section describes our theoretical frame- 
work and discusses experimental limits on the parameters of the model. The following one 
briefly sketches the way we have performed the calculation. Sect. ^ presents some numer- 
ical results, whereas Sect. H summarises our analysis and draws the conclusions. Finally, 
in the two Appendices, we introduce our notation and explicitly derive the Feynman rules 
and cross section formulae needed for our numerical analysis. 
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2 The theoretical model and its parameters 



We work in the theoretical framework provided by the MSSM, the latter including ex- 
plicitly the CP-violating phases and assuming universality of the soft gaugino masses at 
the GUT scale and universality of the soft trilinear couplings at the EW scale. We define 
its parameters at the EW scale, without making any assumptions about the structure 
of the SUSY breaking dynamics at the Planck scale, whether driven by Supergravity 
(SUGRA), gauge mediated (GMSB) or proceeding via other (yet unknown) mechanisms. 
We treat the MSSM as a low-energy effective theory, and input all parameters needed for 
our analysis independently from each other. However, we require these to be consistent 
with current experimental bounds. In fact, given the dramatic impact that the latter can 
have on the viability at the Tevatron and/or the LHC of the CP- violating effects in the 
processes we are dealing with, we specifically devote the two following Subsections to dis- 
cuss all available experimental constraints. The first focuses on collider data, from LEP 
and Tevatron; the second on the measurements of the fermionic EDMs. (Some bounds 
can also be derived from the requirement of positive definiteness of the squark masses 
squared.) Following this discussion, we will establish the currently allowed ranges for the 
Higgs and sparticle masses and couplings. 

Before proceeding in this respect though, we declare the numerical values adopted 
for those MSSM parameters that are in common with the Standard Model (SM). For 
the top and bottom masses entering the SM-like fermionic loops of our process, we have 
used rrit = 175 GeV and = 4.9 GeV, respectively. As for the gauge couplings, the 
strong, electromagnetic (EM) coupling constants and the sine squared of the Weinberg 
angle, we have adopted the following values: as{Mz) = 0.119, aEM^Mz) = 1/127.9 and 
sin^ 9 w{Mz) = 0.2315, respectively. 



2.1 Limits from colliders 

The Higgs bosons and sparticles of the MSSM that enter the gg <l>° production processes 
can also be produced via other channels, both as real and virtual objects. From their 
search at past and present colliders, several limits on their masses and couplings have 
been drawn. As for the neutral Higgs bosons of the MSSM, the most stringent bounds 
come from LEP. For both M^o and M^o these are set - after the 1998 LEP runs - at 
around 80 GeV by all Collaborations p2|, for tan/3 > 10. The tightest experimental 



limits on the squark masses come from direct searches at the Tevatron. Concerning the ti 
mass, for the upper value of tan/? that we will be using here, i.e., 10, the limit on m^^ can 
safely be drawn at 120 GeV or so fairly independently of the SUSY model assumed. 



As for the lightest sbottom mass, m^^, this is excluded for somewhat lower values, see 



^See also [g3| for a more recent and somewhat higher hmit - that we adopt here - on M/jO from 
ALEPH, of about 85.2 GeV for tan/3 > 1 at 95% confidence level (CL), usmg data collected at collider 
centre-of-mass (CM) energies in the range 192 GeV ^ v^ee ~ 196 GeV and a total luminosity of about 
100 pb-i. 
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25[| . Besides, D0 also contradicts all models with < 250 GeV for tan/5 ^ 2, 



Ref. 

A = and fi < (in scenarios with equal squark and gluino masses the limit goes up 



< 260 GeV). 



2.2 Limits from the EDMs 

These are possibly the most stringent experimental constraints available at present on the 
size of the CP-violating phases. The name itself owns much to the consequences induced 
in the QED sector. In fact, to introduce a complex part into the soft SUSY breaking 
parameters of the MSSM corresponds to 'explicitly' violating CP-invariance in the matrix 
elements (MEs) involving the EM current, as the phases lead to non-zero TP form factors, 
which in turn contribute to the fermionic EDMs. In contrast, within the SM, it is well 
known that contributions to the EDMs arise only from higher-order CP- violating effects in 
the quark sector, and they are much smaller than the current experimental upper bounds. 



At 90% CL, those on the electron, |27], and neutron, dn P8[, read as 



|4|exp < 4.3 X 10"^^ ecm , |rf„|exp < 6.3 x lO'^^ e cm . (2) 

As mentioned in the Introduction, if cancellations take place among the SUSY contribu- 
tions to the electron and neutron EDMs, so that their value in the MSSM is well below 
the above limits, i.e., |4|mssm < |4|exp and |(in|MSSM < \dn\e^p, then 0^ and (pA can 
be large. To search for those combinations of soft sparticle masses and couplings that 
guarantee vanishing SUSY contributions to the EDMs for each possible choice of the CP- 
violating phases, we have scanned over the (0^, 0a) plane and made use of the program 
of Ref. 0]. This returns those minimum values of the modulus of the common trilinear 
coupling, |y4|, above which the cancellations work. For instance, in the case of the neutron 
EDMs, the dominant chargino and gluino contributions appear with opposite sign over a 
large portion of the MSSM parameter space. Thus, for a given \fi\, a chargino diagram 
can cancel a gluino one and this occurs for certain values of the gaugino/squark masses 
and a specific choice of 1^41 §]. In general, internal cancellations are more likely among 
the SUSY contributions to the neutron EDMs, than they are in the case of the electron. 
So much so that, in the former case, it is even possible to remain consistently above the 
experimental limits (H) if one only assumes the phase of to be non-zero [Q. 

However, not all the surviving combinations of 0^, 0^ and \A\ are necessarily allowed. 
In fact, one should recall that physical parameters of the MSSM depend upon these three 
inputs. In particular, the squark masses (entering the triangle loops of the production 
processes considered here) are strongly related to 0^, (pA and \ A\. Given the assumptions 
already made on the trilinear couplings (i.e., their universality), and further setting (see 
Appendix A for the notation) 

M,-3 ^ = = , (3) 
M,., = Mq^^^ = M^^^^ = M^^^^ , (4) 
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with Mg^ 2 ~ ^93) where M^^ 2 3 are the soft squark masses of the three generations, one 
gets for the hghtest stop and sbottom masses the following relations: 



nir 



m 



bi 



Ml+ml + -Ml cos2f3 



;^M| - ^M2,)2 cos2 2/5 + 4m2 |A|2 + |/i|2 cot^ /3 + 2|A| cos(0^ - 0^) cot /3 



(5) 



Ml+ml 



-Ml cos 2(3 



■-Ml 
■6 ^ 



-M2 
3 



Wj 



cos2 2(3 + Ami W + lywP tan2 (3 + 2| A| cos( 



tan (3 



(6) 



For some choices of 0^, 0^ and |y4| and a given value of Mgg, and tan/3, the two above 
masses (squared) can become negative. This leads to a breaking of the SU(3) symmetry, 
that is, to the appearance of colour and charge breaking minima. In order to avoid this, 
some points on the plane (0^, 0^) will further be excluded in our study. 



3 Numerical calculation 

We have calculated the Higgs production rates in presence of the CP-violating phases 
exactly at the leading order (LO) and compared them to the yield of the ordinary MSSM 
(that is, 'phaseless') at the same accuracy. In our simulations, we have included only the 
t-, 6-, ti", ^2", ^1- and 62-loops, indeed the dominant terms, because of the Yukawa type 
couplings involved. In order to do so, we had to compute from scratch all the relevant 
analytical formulae for production. In fact, one should notice that for such a Higgs 
state there exist no tree-level couplings with identical squarks if 0/^ = 0a = 0, whereas 
they appear at lowest order whenever one of these two parameters is non-zero. Besides, 
being the quark loop contributions antisymmetric (recall that is a pseudoscalar state) 
and the squark loop ones symmetric, no interference effects can take place between the 
SM- and the SUSY-like terms in the ME for gg A^ . (That is to say that 0^- and 0^- 
induced corrections are always positive if = A°.) The full amplitude is given explicitly 
in Appendix B. For completeness, we have also recomputed the well known expressions for 
scalar Higgs production, <l>° = finding perfect agreement with those already given 

in literature (again, see Appendix B). Here, CP- violating effects can produce corrections 
of both signs. 

It is well known that next-to-leading order (NLO) corrections to gg —>■ processes 
from ordinary QCD are very large [^, |3^. However, it has been shown that they affect the 
quark and squark contributions very similarly [Q. Thus, as a preliminary exercise, one 
can look at the LO rates only in order to estimate the effects induced by the CP-violating 
phases. In contrast, for more phenomenological analyses, one ought to incorporate these 
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QCD effects. We have eventually done so by resorting to the analytical expressions for the 
heavy (s)quark limit given in Ref. These are expected to be a very good approxima- 



tion for Higgs masses below the quark-quark and squark-squark thresholds. Therefore, we 
will confine ourselves to combinations of masses which respect such kinematic condition. 

As Parton Distribution Functions (PDFs), we have used the fits MRS98-LO(05A) §^ 
and MRS98-NLO(ET08) in correspondence of our one- and two- loop simulations, 
respectively. Consistently, we have adopted the one- and two-loop expansion for the 
strong coupling constant as{Q), with all relevant (s)particle thresholds onset within the 
MSSM as described in [Q. The running of the latter, as well as the evolution of the 
PDFs, was always described in terms of the factorisation scale Q = Qp, which was set 
equal to the produced Higgs mass, M$o. In fact, the same value was adopted for the 
renormalisation scale Q = Qr entering the Higgs production processes, see eqs. (|B.1|) and 
(Q-(^). 

Finally, for CM energy of the LHC, we have assumed ^/Sp^ = 14 TeV; whereas for the 
Tevatron, we have taken y/s^p = 2 TeV. 



4 Results 

A sample of 0^, (pA and |y4| values that guarantee the mentioned cancellations can be found 
in Fig. ^, for the two representative choices of Mg^ ^ given in Tab. I. Here, both \^\ and Mq^ 
are held constant at, e.g., 600 and 300 GeV, respectively. The soft gluino mass is given 
too, in Tab. I, as it enters our analysis indirectly, through the EDM constraints (recall the 
discussion in Subsect. |2]^). The allowed values for the modulus of the (common) trilinear 
couplings \A\ are displayed in the form of a contour plot over the (0^,0^) plane, where 
both phases are varied from to vr. (Same results are obtained in the interval (7r,27r), 
because of the periodic form of the SUSY couplings and mixing angles: see Appendix 
A.) In the same plots, we have superimposed those regions (to be excluded from further 
consideration) over which the observable MSSM parameters assume values that are either 
forbidden by collider limits (dots for the lightest stop mass and squares for the lightest 
Higgs mass: see Figs. below) or for which the squark masses squared become negative 
(crosses), for a given combination of the other soft SUSY breaking parameters. Typically, 
we obtain that for small yU phases, i.e., 0^ < 7r/300 ~ 0.01, the value of |y4| tends to 
be zero for almost all values of 0^. In the region where both phases are quite large 
(0At, 0yi ~ ^/2), the modulus of the trilinear coupling must be around 700 GeV, for the 
EDM constraints to be satisfied. 

In Tab. I, in order to completely define our model for the calculation of the gg $° 
processes, we also have introduced the Higgs sector parameters: the mass of one physical 
states, e.g., My^o, and the ratio of the VEVs of the two doublet fields, i.e., tan/3. We 
have fixed the former to be 200 GeV, whereas two possible choices of the latter have 
been adopted, 2.7 and 10 (corresponding to Mg^ j = 1000 and 300 GeV, respectively). 
With regard to this last parameter, tan/5, a few considerations have to be made at this 
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point: namely, concerning the so-called 'Barr-Zee type graphs' Very recently, the 



corresponding contributions to the electron and neutron EDMs have been calculated at 



two loops, in Ref. ^Sf. These terms put bounds (diamond symbols) directly on the 
squark masses and soft trilinear couplings of the third generation and thus are crucial for 
our analysis. However, they are quantitatively significant only at large values of tan /3, as 
can be appreciated in the right-hand plot of Fig. ^. The smaller tan (3 the less relevant 
they are: see the left-hand plot of the same figure. Besides, these two-loop terms entering 
the EDMs also depend on |/x|: small values of the latter induce negligible contributions 
to both the electron and neutron EDMs. (For example, for the choice = 500 GeV 
and tan/9 = 3 made in Ref. [^, one obtains no bounds from the EDMs through the 
Barr-Zee type graphs.) Moreover, for tan/3 ^ 10, in order to be consistent with the 
EDM constraints in (^, one would need the modulus of the soft trilinear coupling to be 
unnaturally large, even greater than 4 TeV (for 0^ ^ cf)^ ^ ^/2). This would drive the 
squark masses in eqs. (|D-(|D to become negative over most of the (0;x,0a) plane. Thus, in 
order to avoid all such effects, we limit our analysis to the interval 2 < tan /5 ^ 10 (recall 
the lower experimental limit on such a parameter: see Subsect. p.l| ). Corresponding values 
for the modulus of the (common) trilinear couplings are in the range |A| ^ 1 TeV. In this 
regime, the parameter combinations given in Tab. I should serve the sole purpose of being 
examples of the rich phenomenology that can be induced by the CP-violating phases in 
the MSSM, rather than benchmark cases. 



tan (3 






M,-3 




Mao 


2.7 


600 


1000 


300 


300 


200 


10 


600 


300 


300 


300 


200 



Table I: Two possible parameters setups of our model. (Apart from the dimensionless tan /3, all 

other quantities are given in GeV.) 



Furthermore, notice that, starting from the numbers in Tab. I, one can verify that 
the heaviest squark masses, m^^ ^'^d m^^, are both consistent with current experimental 
bounds. As for the lightest stop, we display in Fig. ^ the values assumed by m^^ over the 
usual (0^, (pA) plane, for the two choices of tan (5 in Tab. I and in correspondence of the 
\A\ values of the previous figure. As a matter of fact, the effect of the phases is quite 
significant on the actual lightest stop mass, see eq. (^. We observe that, for 0^ ^ |, 
one can get small values for m^-^, if not negative. This is due to the fact that, in such a 
region, \A\ can still get large enough (despite of a low tan/?), so that the last term on the 
right-hand side of eq. (|^) becomes comparable to the first two terms. Anyhow, over most 
of the (0^, 0a) plane, m^^ is well above the current experimental reach of 120 GeV. The 
lightest sbottom mass, m^^, see eq. (H), is always around 290 GeV, so that we have avoided 
to reproduce here the corresponding plots. Also, for the above choices of tan /5, one gets 
constraints on the mass of the lightest Higgs boson. In order to derive these, we make use 
of the two- loop analytic formula for M^o [^. We display the corresponding values over 
the (0^, (Pa) plane in Fig. ^ The mass regions excluded by LEP (see Subsect. |2rT|) amount 
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to a restricted part of the (0^, (pA) plane, both at small and, particularly, large tan (3. As 
for the heaviest Higgs boson masses, one has M^o = 212(201) GeV for tan/? = 2.7(10), 
thus almost degenerate with M^o. 

Having defined the allowed spectrum for the (s)particles masses and couplings entering 
the Higgs production modes that we are considering, we are now ready to quantify the 
effects of the phases on the actual cross sections. A convenient way of doing so is to 
simply look at the ratio between the MEs computed with and without phases. In fact, 
at LO accuracy and further assuming that the relevant hard scale is the same in both 
cases (e.g., Q = M$o), such a ratio coincides with that obtained at cross section level, 
independently of the choice of the PDFs and of y/s. Thus, we define 

- *°) - :ssM ^ . m 

By means of the notation MSSM*, we refer here and in the following to the case of the 
MSSM with same 1^41 but finite values of either 0^ or 0^. Thus, e.g., when 0^ = 0^ = 
(and l^l = too, see Fig. the expression in (0) is of course equal to 1. (Several results 
for this ratio have already been presented in Ref. pO|, though for a choice of the MSSM 
parameters different from those considered here.) 

In Fig. H we present R{gg /i°) as a contour plot over the (0^, 0^) plane for the two 
choices of MSSM parameters of Tab. I. As a consequence of the expression for the gg 
amplitude, see eq. ( p.l|) , the corrections induced by the presence of finite values of 0^ 
and/or (pA in the squark loops can be either positive or negative. Interestingly enough, 
destructive interferences take place over regions already excluded by direct Higgs searches. 
Only constructive interferences would then be observable. These can considerably enhance 
the value of the cross sections obtained in the ordinary MSSM, particularly at large 0a 
and intermediate 0^ values, if tan/5 is small (by up to a factor of 7). This result is 
a consequence of two related aspects. Firstly, in those (0^,0^) regions, \A\ achieves 
its maximum, so that the corresponding Higgs-squark-squark couplings are enhanced 
significantly, with respect to the strength of the Higgs-quark-quark vertices (see Appendix 
A). Secondly, large values of \A\ correspond to small values of m^^ (compare Fig. |^ to 
Fig. ^), this yielding a further kinematic enhancement in the squark loops. For large 
tan/3, the portion of the (0^,0a) plane surviving the EDM constraints is much smaller, 
but the CP- violating effects are still large. For example, in the region Stt/S < 0^ ^ 37r/4 
and 7r/3 ^ 0A ^ 371/4, the MSSM* cross sections can be larger by about a factor of 5 with 
respect to the MSSM ones. 

In Fig. [7| we present similar rates for the heaviest scalar Higgs boson, i.e., $° = 
. Contrary to the previous case, here, one typically obtains a suppression of the 
MSSM* cross sections relatively to the MSSM ones, over the allowed regions of the (0^, 0^) 
plane, both at small and large tan (3. The destructive effect is the consequence of the 
interplay between the Higgs mixing angle a and the minus sign in the third term of the 
relevant XHH^ii coupling, see eq. ( |A.49| ), as opposed to a plus sign in A^jO^j. of eq. ( |A.53| ). 
Quantitatively, the effects of the phases are most conspicuous at small tan (3, when 0^ ~ 
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37r/5 for small (pA values (about one order of magnitude difference between the MSSM* and 
the MSSM, at the most). At large tan/3, over the much smaller (c/)^, 0^) plane surviving 
the experimental constraints, the suppression is at most a factor of 2, when 0^ ~ 47r/5. 

The ratio R{gg A^) is plotted in Fig. As already discussed in Sect. H, here the 



corrections induced by finite values of 0^ and/or 0^ are always positive: see eq. (|B.l 



In fact, the ratio R{gg A^) can become as large as 2 at small tan/?. In this case, the 
maximum is obtained very close to the excluded regions: i.e., 0^ ~ 7r/2 and 0^ ~ 0, vr. 
This pattern can easily be understood by looking at the A^ot^f. coupling of eq. (|A.48|) 
and further recalling eq. ( |A.20| ). (Also notice that, for 0^ = 0^ ^ ^ and tan/5 = 2.7, 
the coupling A^ot^^. becomes zero, hence R{gg — > A^) = 1.) For larger values of tan/5, 
the effects of the phases on pseudoscalar Higgs boson production can even be larger, as 
they are induced through the A^bib^ vertex of eq. ( |A.48| ), which benefits from the tan/5 
enhancement. For example, when 0^ ~ Sn/A and for 0a slightly larger than 7r/2, the 
increase can amount to a factor of 3. 

As general remark on the behaviour of the three ratios, one should notice that they 
are close to unity (i.e., no effects from the CP- violating phases) when 0^ is small for every 
value of (j)A- This can easily be interpreted by looking at Fig. H, since, when 0^ 0, one 
has that |A| — >■ and 0^^; — *■ too (these are the mixing angles of the third generation, see 
Appendix A). Thus, the strength of all Higgs- squark-squark vertices becomes very small 
compared to that of the Higgs- quark- quark ones. The opposite {\A\ -^0 when 0^ — * 
for any 0^) is not true, since here the value is fixed and thus 0^-^ are always non-zero. 

As already intimated in Sect. |^, in order to give realistic predictions for CP- violating 
effects in gg $° processes, one ought to include two-loop QCD effects. We do so 
in the reminder of this Section, by considering the NLO production rates of all Higgs 
states at the mentioned CERN and FNAL colliders, cr^lo^'igg $°). We convert 
the total production cross sections to picobarns and again adopt the input parameters 
of Tab. I. (Incidentally, notice that, being m^-^ > 120 GeV, the relation M$o < 2mq is 
always satisfied, for any g = ti, ^2, ^i, recall the discussion in Sect. |^.) The LHC rates 
are displayed through Figs. p|-p^ whereas the Tevatron ones appear in Figs. |12|-0. We 
present these figures mainly as a reference for experimental analyses. In fact, as far as 
CP-violating effects are concerned, the two-loop QCD dynamics is very similar to the one- 
loop one already discussed. In particular, we have verified that the QCD ii"-factors for the 



quark and squark loops are very similar |3D| over the portions of the MSSM parameter 
space considered here. Thus, to obtain the effects of 0^ and (pA on Higgs production 
via gluon-gluon fusion at NLO, it suffices to refer to Figs. ^|-^ with the normalisation 
of Figs. |9|-|TT| and Figs. |T^-|1^, for the LHC and the Tevatron, respectively. Concerning 
the possibility of actually observing the CP- violating effects at either collider, this is very 
much dependent upon their luminosities. 

At the LHC, with an annual value between 10 and 100 fb~^, all the available (0^, 0^) 
areas can in principle be covered, at both large and small tan /? and for all Higgs states, 
as the production cross sections are never smaller than 5 pb or so, see Figs. p|-pn]. In 
particular, this is true where the effects of 0^ and 0a are larger: compare the areas 
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with high R{gg $ ) values in Figs. ^-§1 to the corresponding rates in Figs. P|-[TT|. For 
example, for $° = h'^, at large (pA, intermediate and for tan/5 = 2.7, cr^lo^* {gg h^) 
is around 200 pb. At large tan/5, when Svr/S ^ (p^ ^ Sn/A and 7r/3 ^ 0a ^ 37r/4, the 
MSSM* cross sections are around 100 pb or more. Similarly, for $° = H^, if tan/? = 2.7, 
when 0^ ~ Svr/S and for small 0^ values, ctJ^lq^ ((7(7 — > if") is of the order of 6 pb. At 
tan/5 = 10, when 0^ ~ 47r/5 (for any 0^), one gets cross section rates around 10 pb. 
Finally, if = at small tan/? and when 0^ ~ 7r/2 with 0^ — 0,7r, one has again 
MSSM* cross sections of the order of 10 pb. For large tan/5, 0^ ~ Stt/A and 0^ ^ ■/r/2, 
(^NLo^*(^^ ^ is about 15 pb. 

Concerning the Tevatron, given the value of integrated luminosity expected at Run 2, 
of the order of 10 fb~^, prospects of detecting CP- violating effects in gg —>■ $° processes 
are very slim. The various cJ^lo^ (99 ~^ $°)'s are notably smaller here, because of the 
reduced gluon content inside the (anti)proton at lower ^/s, for a given M$o value. In 
fact, the production rates over not yet excluded (0^, 0a) regions are never larger than a 
handful of picobarns. One can possibly aim at disentangling CP-violating effects in the 
case of the lightest Higgs boson, at both tan/3's, but only in the usual corners where both 
the corrections and the absolute production rates are largest: see Fig. |12[ As for the 
other two Higgs states, the chances are extremely poor, as cr^lo^'{gg — * is always 
well below the picobarn level, see Figs. p!3Hl^. 

As a final remark of this Section, we would like to mention the following. A peculiar 
feature concerning Figs. |9|-[10|, as compared to Figs, ^j-^, respectively, is the different 
pattern of the level curves. This should not be surprising though, as, for = h'^, H^, the 
squark loop contributions in the ordinary MSSM are non-zero to start with. In contrast, 
one can appreciate the strong correlations between the level curves in Fig. 11 and Fig. |^, 
a consequence of the absence of scalar loops in gg ii (pn = (pA = ^- Similarly, for 



the case of Figs. 12-14 



5 Summary and conclusions 

It is well known that finite values of the mixing angles d^j^, converting the weak into the 
mass basis of the third generation of squarks, see eqs. ( |A.6| ) and ( [A.9|) , imply that left-right 



chiral currents, eqs. ( [A.26D -( [A.37[ ), can enter the Higgs-squark-squark couplings appear- 
ing in the scalar loops contributing to gg — *■ processes, if $° = /i^, H^, because of the 
structure of the mixing equations ( [A.ISI ). As a consequence, the corresponding $° pro- 
duction cross sections develop a dependence on /i and A, the Higgsino mass term and the 
trilinear scalar coupling (the latter assumed here to be universal to all (s) quark flavours) 
entering the soft SUSY breaking sector. The strength of their contribution is however 
modified if CP-violating effects are manifestly inserted into the MSSM Lagrangian, by 
allowing these two parameters to be complex, see eqs. ( |A.48| )-( [A.56| ). In such a case, in 



particular, also the gg A^ cross section receives scalar loop contributions, thus acquir- 
ing a dependence upon fi and A, much on the same footing as when $° = h^, if": see 
eq. 
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Clearly, it is the actual size of the independent phases associated to the above two 
parameters, (p^ and 0^, that regulates the phenomenological impact of complex values of 
fi and/or A, not least, because they also affect the two mixing angles: see eqs. ( |A.19|) and 
( A.24| ). Given the importance of Higgs production at future hadron-hadron colliders, such 
as the Tevatron (Run 2) and the LHC, and the fact that gluon-gluon fusion is a sizable 
production mode at the former and indeed the dominant one at the latter (over most of 
the MSSM parameter space), we have made the investigation of the effects induced by 
finite CP-violating phases in gg $° processes the concern of this paper. In order to 
address the problem quantitatively, we first had to derive the relevant Feynman rules of 
the MSSM in presence of (j)f^ and (j)A (Appendix A) and eventually calculate the associated 
cross sections, for any Higgs state (Appendix B). 

Before proceeding to the numerical analysis though, we had to introduce a parametri- 
sation of our theoretical model and incorporate the latest experimental constraints on its 
parameters. These can be subdivided into two categories, those arising from analyses per- 
formed with collider data and those deduced from the measurements of the electron and 
neutron EDMs. The former mainly limit the value of the squark masses and couplings, 
thus only indirectly affecting 0^ and (j)A, see eqs. (il)-(|). In contrast, the latter can be 
very stringent in this respect, unless cancellations take place among the SUSY contribu- 
tions to the fermionic EDMs, so that the bounds derived this way on 0^ and 0a become 
much less potent. In practise, the CP- violating phases can attain any value between 
and TT, provided and |y4| are in appropriate relations, which are in fact satisfied over 
large portions of the MSSM parameter space. Under these circumstances then, the CP- 
violating phases can affect the interplay between the quark and squark loops in gg — »• 
processes considerably (Section 2). 

In the end, we have verified, both at LO and NLO accuracy (Section 3), that this 
is true, over those parts of MSSM parameter space where these cancellations are more 
effective. As a matter of fact, effects due to finite values of 0^ and/or 0^ can be extremely 
large, inducing variations on the Higgs cross sections of the ordinary MSSM (i.e., those 
obtained for 0^ = 0yi = at the same \A\) of several hundred percent, at least for 
values of tan/? in the range between 2 and 10 and soft masses and couplings below the 
TeV region (Section 4). For these combinations of parameters, even the bounds induced 
by the contributions of the Barr-Zee type diagrams to the EDMs can easily be evaded. 
Other than studying relative effects of the phases, with respect to the yield of the ordinary 
MSSM, we also have presented absolute rates for the gg — >• $° cross sections at NLO, 
at both the LHC and the Tevatron, for the MSSM including 0^ and 0a, thus showing 
CP-violating effects explicitly in observable quantities. Given the higher luminosity and 
production rates expected at the CERN collider, as compared to the values at the FNAL 
one, real prospects of sizing these effects in Higgs boson production will most likely have 
to wait for a few more years. By then, the available portions of the (0^, 0^) plane should 
also be expected to be better defined than at present, given the improvement foreseen in 
the near future in the precision of the EDM measurements [37|. 

Anyhow, as we have tried to motivate in the Introduction, and following our results. 



11 



we believe that further investigation is needed of the consequences of exphcit CP-violation 
being present in the soft SUSY Lagrangian. For example, to stay with the Higgs sector, 
one should establish the effects of (p^ and (pA in the decay process hP — > 77 ||39|, as this 



represents the most promising discovery channel of the lightest Higgs boson of the MSSM 
at hadron-hadron machines. In this case, the proliferation of SUSY induced contributions 
(also due to charged Higgs bosons, sleptons and gauginos) could well be responsible of 
CP- violating effects comparable to those seen in the Higgs production processes via gluon- 
gluon fusion, as the latter are solely due to squarks contributions. 
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Appendix A: the CP-violating phases in the MSSM 

In this Section, we follow the notation of Ref. ||3^. We start from the Superpotential, 
which has the form 

W = e,j{YeHiL^E + YdHlQW + Y^HiQ'U + fxHlH^) (A.l) 

and where all fields appearing are actually Superfields, with ei2 = 1. In terms of compo- 
nent fields, the Lagrangian of the soft breaking terms reads as 

Csoft = -Ml^ liJip - M^JifaP - ^iBei, {H{Hi + h.c.) 

- M? (uluL + dldL) - M^u^ur - Mld^Dn 

- e,, (-YuA^HlQ^u*R + Y^A^HlQ^d^R + YAeH{Ue*R + h.c.) . (A.2) 
The squark mass squared matrix (here and in the following, g^'^ = t and h) 

is Hermitian and can be diagonalised by the unitary transformation 

UlMlU, = diag(M|^,M|J, (A.4) 

with (M,-2 < M^2) 

Km = lii^L + MIrr) (-)[+]v/(M|,, - Ml^^y + 4 \M,n,\^] (A.5) 



and 



cos Or, — sin Or, e**^? 
sm 6q e cos 6q 



where — 7r/2 < 9q < rt /2 and 



2 |M?R.| 

tan = ' , (A.7) 



MIl - M?RR 



2 
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where refers to the imaginary part of a complex quantity. 

Now, in order to construct the Feynman rules for the Higgs-squark-squark vertices, 
we need the transformation of the squark weak basis {qi, qn) into the mass basis (gi, ^2)5 
namely 

''')^uJ':). (A.9) 
qn J J 

Furthermore, one has to proceed by also transforming the Higgs boson weak ba- 
sis [H^, H2, H2 , H2) of the four complex fields into the real eight physical ones 
{H°,h^,A^,H+,H-,G^,G+,G-). Following Ref. the transformation can be writ- 
ten as follows 

= t;i + -^(i7°cosa-/i°sina + iA°sin/?-zG°cos/?) , (A. 10) 

= W2 + ^ sin a + cos a + cos /5 + sin /3) , (A.ll) 
v2 

H{ = H- sin f3 - G- cos (3 , (A.12) 
H+ = H+ cos (3 + G+ sin (3 , (A. 13) 

with {H+y = H" and 



where M^o , M/jO are the tree- level CP-even Higgs boson masses. 

The Feynman rules for the Higgs-squark-squark vertices, involving mixing and phases, 
finally are (here, $ = $° and H^): 



"^'I'g2gi* !^qW^ '^^qLq£ ^ ^q^q'^ '^^qm^ i'q^'q'fi '^^qL'l^ ^ ^<S>qaq'£ ■ 



(A.15) 



For the case of stop squarks (i.e., q = t), one has 



ML = M?+M2 + -(4M2,-M|)cos2/3, (A.16) 
MLr = M? +M2-^(M^-M|)cos2/3, (A.17) 



3 



MinL = iMiLRr = mt{A + fi*cotf3), (A.18) 
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and thus 



t3Jl26z = 5 r-T , A. 19 

M|-M? + (|M2,-|M|)cos2/3' ^ ^ 



sin^t = ^j— — - — —T, , (A.20) 





sin 


- 


l/^l 


sin 0^ cot /? 




\At + fi* cot/3| 



where fx = |/i|e*''^^ and = \At\e^^^-t. 

Similarly, one obtains for sbottoms (i.e., q = h): 

^III = M?+m2-l(2M^ + M|)cos2/5, (A.21) 

MLr = M? +m^ + i(M^-M|)cos2/5, (A.22) 

Ml^^ = {M?^^y = mh{Ah + fi*ta.n(3), (A.23) 

with 

tan26'r = — — — — rr , A. 24) 

f M2-M2 + (-|M2, + iM|)cos2/5' ^ ' 

. |y4fe| sin0^^ - sin0^tan/3 /aoc;\ 
\Ab + fx* tan PI 

where Ai, = \Ab\e''^'^b. 

The chiral couphngs A^^^^', (x = L,R) of eq. ( |A.15| ) can be found in Ref. with 
the only exception of those cases where fi and Ag enter^, for which one has to adopt the 
following set of formulae (with g'^ = Ana-Eu/ sin^ 9w the weak constant and where Mw is 
the boson mass): 

Aao*,*- = (/i* - A^ cot P) , (A.26) 

^AOiliu = -i^AOiLi}J* ) (A.27) 

■^^0^*^^ = -{^Ao'bLb'J* ' (A.29) 

^HHltn = -i^H"iLi*J* 5 (A.31) 
•^In other terms, for real /i and Aq our expressions reduce to those in Ref. ||3^ . 
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^H%lhR — ~i^HObLb*J* ' (A.33) 

= 2M^Li5 " " ' ^^-^^^ 

^hOiliu = -{^hOiLi'J* ) (A.35) 

Koi^b* = -7TTr^^Z^if^*co^(^~'^dSina) , (A.36) 

^ ZlVlw COS p 

'^hOblbn = 'i-^hO'bLb'J* ' (A.37) 

= - /3) , (A.38) 

= -^|^(/i-A*tan/3) , (A.40) 

(Although we have not made use of the Feynman rules involving charged Higgses we 

display them here for completeness.) We need also the interactions among Higgs bosons 
and quarks and these read as follows: 

qrriu cot Q , , . 

^AOif = - ^ 75, (A.42) 



2M, 



w 



^AObb = ^ 75 , (A.43) 



w 



2Mw sm p 

igrudcosa 
2Mw cos p 

igmucosa 
2Mvy sm p 

_ _ ^^mrfsino; 

2Mw COS p 

It is now useful to look at the explicit phase dependence of the vertices involving Higgs 
bosons and squarks: to this end, we expand our formulae ( |A.15| ). The relevant couplings 
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involving the CP-odd Higgs boson read as^: 

-n2%||/x|sin((/)£-(/)/,) - |^t|sin((/)f + (/)Ajcot/3| 



.2Mw- 



sm{4>i - (j)^) - \Ab\ sin 



For the CP-even Higgs bosons we find (here, sw = sin 9w), 



#6 + <PA, 



) tan/Sj , 



(A.48) 



"^HO'bibl 



"^ROhb^ 



'igMz' 


{-[ 


Cw 





1 777-^ sin cx 

- cos^ 9i - euSw COS 2% cos(a + p) - -jA^^ 
1 -I M| sm p 



rrit sin 2% 
2M| sin /3 





{-[ 







cos(0f — (j^^) cosa + |A| cos{(f)i + ^^t) sin a 

9 / /i\ ""^^ sin a 

cos 26^ cos a + /?) - , ■ a 
J M| sm p 



, (A.49) 



^ sin^ 0f + e„s^ cos 



sin 2^1 r, . , , , , , ^ . , , , s . 1 1 / » ,-^\ 
i"i cos(0f — cosa + |74t| cos(0f + 0aJ smo; > , (A. 50) 



rrit sin 2^^ 
^ 2M|sin/3 



ic/M^l jrl ^ 2 o/j1 



nib 



^, m? coscK 
^ Ml cos 13 



'tgMz' 


{[ 


Cw 





1.2/1 2 ^/i 1 / /,\ cosa 

- sm - edS^^ cos 2^ cos(o; + /3) - ^^^^ 

sin 26',, r, I , , , , . , . , 
+ cos/3 r ' ~ '''^ " ^ ' ' 



igMz rr- 



x 

- cos^ - e„s^ 
Cw Lz 

sin2^,-r, , 



•s{(f)i + 0aJ cos a 
cos a 



, (A.52) 



mtsin26'f r 
^ 2M|sin/5 



'igMz' 


{[ 


Cw 





cos(0j — 0^) sin q; 
- sin^ 9^ + e„s^ 



2 0/1 1 • / "^f cosa 

-2«,-Jsm(a + /3)-^— 

— \At\ cos{4>i + 4>At) COS a 



, (A.53) 



cos 29^ sm{a 



, cos a 



''Note that the vertices with gluons and squarks are not affected by the phases. 
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rrit sin 26 1 
2M|sin/3 L 



cos((/)(- — (/)^) sin a — l^tl cos((/)(- + (/)^J cos a >, (A. 54) 



A 



h%ibl 



'igMz' 


{-[ 


Cw 





ml sin a 



- cos^ Oi + CrfS^ cos 20^J sin (a + ' ^2 ^ 



mfe sin 26*^ 
2M| cos (3 



cos(0^ — (j)^) COS a — l^bl cos(0^ + (pA^) sin a > , (A. 55) 



A 



•202 



'igMz' 


{-[ 


Cw 





+ 



mf, sin 26*^ 
2M| cos /5 



X 1 t/1l sill 

- sin^ $1 - e^s^ cos 2^^J sin(a + 13) + 

li\ cos(0^ — 0^) COS a — \ A]\ cos(0^ + sin a 



(A.56) 



where e„ = +2/3 and = —1/3. Of course, one can get the corresponding vertices for 
all squarks flavours by a simple substitution t — >■ m, c and b ^ d,s. 



Appendix B: matrix elements and cross sections 

Since, as we have already discussed in the main body of the paper, the CP-violating phases 
induce a non-zero contribution from squark loops in CP-odd Higgs boson production which 
is absent in the phaseless MSSM and since this has not been evaluated yet in the literature, 
we had to perform the loop tensor reduction in such a case from scratch. However, for 
comparison purposes, we have also recalculated the well known tensor associated to CP- 
even Higgs boson production and found agreement with old results. We have evaluated the 
diagrams in the (modified) Dimensional Regularisation {DR) scheme, which preserves the 
SUSY Ward identities of the theory up to two loops. This enabled us to check analytically 
the gauge invariance of our results. Furthermore, we have carried out the 7-algebra in 
four dimensions while using analytical continuation in d dimensions in order to calculate 
the divergent parts of the integrals. The squared matrix elements summed/averaged over 
final/initial spins and colours finally areQ: 



2567r2 



E 



A 



hOqq 



1 + 



1 A/jOgq* 



1 - r,-/(r,-) 



HO 



2567r2 
2567r2 



rua 



m 

^AOqq 



1 + (1 - r,)/(r,) 



1 ^HOqq* 



1 ' 1 



TqfiTq 



^lelE 



m- 



B.l) 



^Note that the 75-matrix is here intended to be removed from the expressions of the A^o^^'s of 
eqs. ( |A.43D -( |01 . 
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4mt 



where Tg^g — -jtt^, q — t,b and q — ti, t2,bi,b2- The function /(r) stands for 





4y(i - y) 



arcsin2(^), 



r > 1 

T < 1 



(B.2) 



and some useful hmits are 



^lim{r[l + (l-r)/(r)]}=+-, 
lnn{.[l-./(.)]} = -l, 



lim s r J 

r— *oo ■ 



(B.3) 
(B.4) 
(B.5) 



The non-existence of interference terms between quark and squark loops for CP-odd Higgs 
boson production can readily be understood by looking at the corresponding amplitude 
formula (Pi, P2 are the gluon four- momenta, e^(Pi), e,y(P2) their polarisation four- vectors 
and a, b their colours) 



as{Q) 
27r 



Sabe^^{Pl)e^{P2) x 



(B.6) 

where one notices an antisymmetric part - note the Levi-Civita tensor e - associated to the 
quark contributions (first term on the right-hand side) and a symmetric one associated 
to the squark loops (second term on the right-hand side). For completeness, we also 
give the tensor structure of the loop amplitudes corresponding to CP-even Higgs boson 
production: 

ie,(Pi)e.(P2)M^'^(g9 - h\ H') = ^<5„,e^(Pi)e.(P2) [g^'^P, ■ P2 - P^,P^) x 



m„ 



1 + (1 - Tg)f{Tg 



27r 
1 

" 4 



E 



A 



{h°,HO)qq* 



1 - rgf{Tg) 



(B.7) 



Here, interference effects clearly exist between the SM- and SUSY-like parts, because of 
the symmetric nature of both contributions. 

The LO partonic cross sections at the energy are then 



s 



TT 
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whereas the corresponding hadronic rates for a colhder CM energy -^i read as 

TT ,-r-T,2 dC 



4o = r ^ , (B.9) 



with r = - = and 

s s 

r^ = l'^rg{x,Q)gC-,Q), (B.IO) 

where g{x, Q) is the PDF of the gluon, evaluated at the scale Q = M$o. 

As already mentioned in the Introduction, in performing the two-loop analysis, we 
have made use of the formulae given in Ref. We refer the reader to that paper for 
specific details. 
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Figure 1: SM-like contributions from top (t) and bottom (b) quarks to Higgs boson 
production via gg $° in the MSSM. 




Figure 2: SUSY-like contributions from top (^1,2) and bottom (61,2) squarks to Higgs boson 
production via gg — > $° in the MSSM. (Notice that, if the CP-symmetry is conserved, 
then $V ^°-) 
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Figure 3: Contour plots for the values of the modulus of the common trilinear coupling, 
l^l, needed in order to obtain the cancellations of the SUSY contributions to the one- 
loop EDMs, over the (0^, 0^) plane for small (left-hand plot) and large (right-hand plot) 
tan/3. The other MSSM parameters are as given in Tab. I. Here and in the following, "x" 
symbols denote points excluded because of the negativity of the squark masses squared; 

symbols denote points excluded by the two-loop Zee-Barr type contributions to the 
EDMs; "□" and "•" symbols denote points excluded from Higgs boson and squark direct 
searches, respectively. 
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ml (GeV), tan/3 = 2.7 



[rad] 
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m^^ (GeV), tan,8=10 
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xxxxxxxxxxx 
■ xxxxxxxxx 
xxxxxxxxxxx - 



[rad] 



Figure 4: Contour plots for the values of the lightest top squark mass, m^^, corresponding 
to those of \A\ in Fig. 0, over the {(p^^cpA) plane for small (left-hand plot) and large 
(right-hand plot) tan j3. The other MSSM parameters are as given in Tab. I. 




Figure 5: Contour plots for the values of the lightest Higgs boson mass, M^o, correspond- 
ing to those of \A\ in Fig. ||, over the (0^, 0a) plane for small (left-hand plot) and large 
(right-hand plot) tan (3. The other MSSM parameters are as given in Tab. I. 
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[rad] 




R(g g ^ hO), taii,S=10 



[rad] 




Figure 6: Contour plots for the values of the ratio in eq. (|^ for the case $° = h^, 
corresponding to those of 1^41 in Fig. over the (0^,0^) plane for small (left-hand plot) 
and large (right-hand plot) tan (3. The other MSSM parameters are as given in Tab. I. 
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R(g g ^ HO), tan|S=10 
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Figure 7: Same as in Fig. ^ for the case $° = if °. 
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Figure 8: Same as in Fig. ^ for the case $° = A^. 
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Figure 9: Contour plots for the values of the NLO cross section for gg $ in the MSSM* 
at the LHC, cf^Iq^' {gg — > $°), for the case $° = corresponding to those of \A\ in 
Fig. 1^, over the (0^, 0^) plane for small (left-hand plot) and large (right-hand plot) tan {3. 
The other MSSM parameters are as given in Tab. I. 
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a(g g ^ HO) (pb), tan,S = 2.7 ^(g g ^ HO) (pb), tan/S=10 
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Figure 10: Same as in Fig. ^for the case $° = H^. 



(j(g g ^ AO) (pb), tan|S = 2.7 ^(g g ^ AO) (pb), taii|S=10 
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0/2 [r^^] ^t^ij^ [rad] 

Figure 11: Same as in Fig. ^for the case $° = . 
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a(g g ^ hO) (pb), tan,8=2.7 cj(g g ^ hO) (pb), tan,S=10 
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[rad] [rad] 



Figure 12: Contour plots for the values of the NLO cross section for gg in the 

MSSM* at the Tevatron, cr^^o^* {gg — ^ $°), for the case $° = corresponding to those 
of \A\ in Fig. over the plane for small (left-hand plot) and large (right-hand 

plot) tan /3. The other MSSM parameters are as given in Tab. I. 



a(g g ^ HO) (pb), tan,8 = 2.7 



a(g g ^ HO) (pb), tan|S=10 



3 —I 



[rad] 



□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□□□ 
□□□□□□□□□□□ 
□□□□□□□□□□□ 
□□□□□□□□□□□□ 
□□□□□□□□□□□□ 

(^□□□□□□□□□n 
□□□□□□□□□□ 
□□□□□□□□□□□ 
□□□□□□□□□□□ 
□□□□□□□□□□ 
□□□□□□□□□ 

□□□□□□□□ 

□□□□□□ 




□ □□□ 
□□□□□□□ 
□□□□□□□ 
□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□ — 
□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□ □□□□□□□□I 



(p., [rad] 




Figure 13: Same as in Fig. Ofor the case = H^. 
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Figure 14: Same as in Fig. 12 for the case $° = A^. 
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